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ABSTRACT 

We study the phenomenology of the Peccei-Quinn invariant extension of the next-to- 
minimal supersymmetric standard model (NMSSM) in view of the recent discovery of a 125 
GeV Higgs boson. The minimal model having no quadratic and cubic terms of the NMSSM 
singlet field predicts a light singlino-like lightest supersymmetric particle (LSP). The model 
is strongly constrained by the Higgs invisible decay and the dark matter characteristic of 
the LSP, while some constraints can be relaxed by assuming that the saxion, the CP-even 
companion of the axion in the Peccei-Quinn sector, causes a late-time entropy production 
diluting the thermal LSP density. The collider signal of the model contains multi-jet and 
h/W/Z plus missing energy, which can be discovered in the early stage of the 14 TeV LHC 
running. 
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1 Introduction 



The strong CP problem can be nicely resolved in a supersymmetric (SUSY) standard model, 
which allows an extended Higgs sector to implement the Peccei-Quinn (PQ) symmetry [1]. 
After integrating out the heavy PQ sector around vpq = 10^ — 10^^ GeV, and freezing 
the axion supermultiplet at their vacuum values, the low energy theory can be reduced to 
the conventional minimal supersymmetric standard model (MSSM) [2], or to the next-to- 
minimal supersymmetric standard model (NMSSM) [3] as proposed recently in [HE]. 

In the PQ-symmetric NMSSM (PQ-NMSSM), the low-energy theory is generically 
described by the superpotential, 

W = XSH^H, + fxlS + ^fx'sS^ (1) 

and corresponding soft-breaking terms in the scalar potential, where fis s-nd fig at the TeV 
scale can arise as a function of vpQ after the PQ symmetry breaking at the scale vpq. 

In this paper, we investigate the phenomenology of the PQ-NMSSM in light of the 
recent LHC results on the Higgs boson search [71[8]. In its minimal form [1], the model 
contains only the /x| term predicting a very light singlino. Because of this, the minimal 
model is severely constrained by the consideration of the Higgs invisible decay and the dark 
matter property of the singlino-like lightest supersymmetric particle (LSP) in the standard 
cosmology. A sizable NMSSM contribution to the 125 GeV Higgs boson mass can be 
obtained with tan /3 ~ 1, which makes the LSP heavy enough to forbid the invisible Higgs 
decay. However, the recent XENONIOO result on the direct detection of dark matter [9] 
excludes almost all the LSP mass region, although the LSP annihilation by the Higgs 
resonance effect can reduce the dark matter relic density significantly. For larger tan/3, 
the LSP gets lighter to open the Higgs invisible decay channel. In this case, the coupling 
between the Higgs boson and the LSP can be made small by a cancellation. Considering 
the charged Higgs boson mass bound from the LHC, the Higgs invisible decay branching 
fraction can be smaller than 0.1 for tan/3 > 9. For such a large tan/3, the NMSSM 
contribution to the Higgs boson mass becomes negligible. As the light singlino being the 
LSP couples very weakly to the quarks and leptons, its thermal relic density overdoses 
the universe. We will argue that the difficulties of the thermal LSP dark matter can be 
circumvented by a late-time entropy production by the saxion field which is inherent in the 
PQ-NMSSM, and still a correct amount of dark matter can be provided by either axion or 
non-thermal LSP. Analyzing the generic collider signatures of multi-jet and h/W/Z plus 
missing energy for two chosen benchmark points with small and large tan /3, we find that 
a 5a discovery can be achieved in the early stage of the 14 TeV LHC running. 

This paper is organized as follows. In Sec. [2], we provide a general description of 
PQ-NMSSM. Then, we analyze the Higgs and the neutralino sectors of the minimal PQ- 
NMSSM. In Sec. [3l we study the phenomenology of a 125 GeV Higgs boson and a light 
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singlino to constrain the model parameter space from invisible decays of the Z and the 
Higgs boson, and the direct production of light neutralinos at LEP II. Analyzing the dark 
matter property of the LSP as a mixture of the singlino and Higgsino, the minimal model 
will be tightly constrained in Sec. HI In Sec. [5l we analyze the collider signatures of the 
PQ-NMSSM at the LHC. We conclude in Sec. El 



2 Peccei-Quinn Symmetric NMSSM 
2.1 Model of PQ symmetry breaking 

The U(l)pQ symmetry, which solves the strong CP problem of the standard model (SM), 
should be spontaneously broken between 10^ and 10^^ GeV. In SUSY models, it can be 
realized by introducing PQ charged but SM singlet chiral superfields 

Xi = 4>i^^f2dlii^d''F''\ (2) 

whose scalar components are stabilized at an intermediate scale, (10/1) ~ upq ~ 10^ — 10^^ 
GeV. Then, the QCD axion corresponds to a linear combination of the axial components 
of 0/. Such a large vacuum value, compared to the weak scale, can be easily obtained if 
(10/ 1) are determined by the interplay of soft SUSY breaking terms and the F-term scalar 
potential suppressed by a cut-off scale. In order to write down the higher dimensional term, 
we introduce two PQ-charged chiral superfields and the superpotential, 

l^PQ = ^XrX,4-" (n = l,2), (3) 
Mpi 

where Mpi ~ 2.4 x 10^*^ GeV is the reduced Planck mass. The PQ symmetry is realized as 
U(l)pQ : Xi e-^^-i^Xi, X2 ^ e-^^-2"X2, (ng^, + (4 - n)qx, = 0). (4) 
When one includes the soft SUSY breaking terms, the scalar potential is 

= <i<^ir+<j02r+(^^^0i0r"+h.c.) 

All the soft parameters l^xj ~ |'^X2l ~ I^pqI of the order of mgoft = C(10^ — 10'^) 
GeV. In the case of n = 2, the additional Z2 symmetry [Xj — )■ {—l)^Xj) is needed to 
prevent MxXiX2 term in the superpotential. The PQ symmetry breaking scale can be 
estimated as 

^'PQ = a/^(4(I0iP) + #(I02P)) ~ (|0i|) ~ (|02|) ~ y/m3oftMpi/«:,^, (6) 
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and it is naturally lying on the axion window for a moderate value of k^^. Non-zero auxiliary 
F-components of Xi and X2 are also developed. In canonical basis of the superfields, they 
are given as 



frisoif (7) 



)i / \(pi u(p. I \ 101 P / Mpi 



"^soft, 



^2 



|2 



After fixing the vacuum values, the scalar fields can be decomposed as = ( | + 
(s/ + iai) Masses of the PQ fields (s/, a/, and a/) are generically of order of mgoft 

except the mass of the QCD axion oqcd = ^f^'^vq, of order y'm^A^^/fpQ, 

where is a light quark mass. The PQ charged particles with masses of (9(msoft) can 
play a important role in cosmology, which will be addressed in Sec. HI 

The axion solution to the strong CP problem can be realized by extending the Higgs 
sector through a superpotential term, XlHuHd/M-p\ with qHuH^ = — 2gxi [2]. A simple 
consequence of this extension is that a bare Higgs-bilinear term is forbidden by the PQ 
symmetry and the right size of the yii-parameter is generated, 

The low energy theory after the PQ symmetry breaking could also be of the NMSSM 
type [Hin] where the /iq term is extended to a scalar field S and the Higgs boson mass at 
125 - 126 GeV could be realized with less fine-tuning [TOl[TT]. In the following subsection, 
we will describe how the PQ symmetry is incorporated into the NMSSM setup. 



2.2 General PQ-NMSSM 

In the NMSSM, the SM singlet superfield S couples to the Higgs fields with a renormalizable 
term XSHuHd in the superpotential. Then, 5* should have the PQ charge, qs = —qn^H^- 
Below the PQ symmetry breaking scale, the effective theory will be described by the general 
NMSSM dU: 

IVeff = (/io + ^S)H^Hd + ^llS + ]^^i'sS^ + ^-t^S\ (9) 

where //q, /i|, /^g, and n are the effective parameters determined by interactions between S 
and Xj. Obviously, they should be vanishing in the PQ symmetric limit, wpq — )■ 0. In order 
for S* to be a light degree of freedom that survives around TeV scale, /io, /is, /i'^ < O(msoft), 
and corresponding soft SUSY breaking parameters are around rrisoit as well. Note that /xq 
and the corresponding soft SUSY-breaking parameter (6'^ component of the spurion field 
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Ho) can be always rotated away by the holomorphic field redefinition 5 — j- 5 — /io/A. So we 
can set /io = without loss of generality. Now it is expected that sizable values of /x| and 
fi'g can be obtained by the PQ symmetry breaking model. We will study explicitly how 
such PQ symmetry breaking parameters can be generated for the models given in Sec. 12.11 

For the PQ symmetry breaking model (|3]) with n = 1 (Wpq = Hp^^XY^ /Mp{), the 
relevant PQ invariant superpotential is 

W = + XS^ H^H, + ■ ■ ■ , (10) 

where dots denote highly suppressed terms that do not contribute to the Higgs phenomenol- 
ogy. The non-trivial /i| can be obtained from a PQ invariant higher dimensional term in 
the Kahler potential, 

K ='^XfS + h.c.. (11) 
Mpi ^ ^ ' 

At low energy, kxs and Xj can be regarded as SUSY breaking spurion superfields with 
Kxs = + e'^Kprnsoft + O'^Kprnsoh + 6'^6'^K,Dm'^oit^ and Xj = (0/)(l + e'^{F^' /(pi)), kq, kf, 
Kp, kd are all 0{1) constants with a reasonable assumption that they and the soft SUSY 
breaking terms for the superpartners of the SM fields have the same origin. Then, by ([6]) 
and (CD, 

' Ut 



but still fi'g and k are suppressed. Again note that the superpotential term XfHuHd/Mpi 
in (fTOj) can be removed by the holomorphic field redefinition S ^ S — Xf/XMp\ without 
loss of generality. Then, the low energy effective superpotential corresponds to the minimal 
type of PQ-NMSSM [12], 

W,s = XSH^Hd + filS. (13) 



As for the model with n = 2 (Wpq = Kp^X^F^/Mpi), the same superpotential ffTOl) 
and Kahler potential ( fTTj) are allowed so that the sizable /iq, /i| are generated. As discussed 
in Sec. 12. 1^ there is an additional Z2 symmetry Xj — t- { — 1)^Xj, S — —S, HuHd — )■ —HuHd 
at renormalizable level, in order to insure that terms like X1X2, X2S are absent. Since the 
Z2 symmetry is explicitly broken by the term XfHuHd/Mpi, a cosmologically dangerous 
domain wall is not produced. Besides, the Z2 breaking tadpole induced by supergravity 
loop corrections are suppressed due to the PQ symmetry. This model is more complicated 
than the model with n = 1. However, a sizable /ig term can be generated if there are matter 
superfields Zi, Z2 whose masses are given by (l^il) from a superpotential term X1Z1Z2 with 
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the PQ charges qz^ = 3gxi, <lZ2 = ~4gxi- Under the Z2 symmetry, Zi — )■ {—l)^Zj. Then, 
the following superpotentials are allowed 

AW = \xzX,Z,Z2 + \zsZ2S' + -T^ UxzX^Z, + ^xzsX'.Z^s) + ^X^S + ■■■ . 

Mpi V / Mpj 

(14) 

Since the masses of Zi and Z2 are of the order of vpq, they should be integrated out, at 
the PQ symmetry breaking scale, by the superfield equations of motion: 

0. (15) 



dZi dZ2 

By substituting the solutions of f|T5l) to Zj of f lT^ . A/i| and /i^ are generated as 
.2 _ / Kxzs^xz{<l)i) ^ ^ \ f (02) ^ \ Y ^'PQ ^ ^ 



This case corresponds to the singlet extension of the MSSM [13] . 

l^efr = XSHuHa + + ^l'sS^ (17) 

while K still suppressed and rotated away. The suppressed k is a generic consequence of 
the PQ extension of the NMSSM. 

The PQ sector contributions to the soft SUSY-breaking parameters for Higgs and 
singlet sector are coming from 6*^ {O"^) component of the spurion superfields Xi = (</>/) (1 + 
9'^{F^^ /(f)j)) (Xj), which are the same order of the MSSM soft SUSY-breaking parameters. 



2.3 Minimal PQ-NMSSM 

In this subsection, we investigate the phenomenological consequences of the minimal PQ- 
NMSSM [1] whose low-energy effective superpotential takes the form of (fT3|) in addition 
to the usual MSSM Yukawa superpotential terms. The new superfield S is singlet under 
the SM gauge group, and it can acquire the vacuum expectation value (VEV) to give the 
natural size of the /i term of the electroweak (EW) scale. 

2.3.1 Higgs sector 

Here, we describe the Higgs scalar potential, its vacuum structure, and the mass spectra 
around the vacuum. The Higgs scalar potential consists of the following F- and D-term 
contributions from the superpotential f|T3l) . and soft SUSY-breaking terms, 

vf = \\KH- - + /.If + + mi' + ii/.T + m') , (is) 
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vo = ^ + - |/7,T - m'T + 1 ih^ht + //x-t ^ (19) 

+ [XA^iH^H- - HlH',)S + tsS + h.c] . (20) 

We assume that all the coefficients in the potentials are real so that any explicit CP 
violation does not occur other than the Cabibbo-Kobayashi-Maskawa phase. The VEVs of 
the charged Higgs field should be vanishing in order to obtain a successful EW symmetry 
breaking minimum. One of the VEVs of the charged Higgs field, e.g., {H^) can be made 
zero with positive = Vu by the SU{2) gauge choice. The other one, however, is not 
guaranteed to vanish in contrast to the MSSM case in which the minimization condition 
dV/dH^ = gives (-f^j) = 0. Still, it can be shown that (H^) = satisfies the extremum 
conditions and it will turn out to be local minimum if all masses squared of the charged 
Higgs sector are non-negative. Therefore, it is reasonable to assume that {H~) = and 
the charged Higgs boson masses squared are required to be positive. The Higgs potential 
can then be written as 

W = + W + « + m - [KAxS + f^DH^H', + h.c.] 

{w - + A' i^Mr + -ii^i' + [tss + h.c.] + 4. (21) 

Expanding the neutral Higgs fields around their VEVs, one gets 

- , (22) 

- -j= , {26} 

Sr + iSi 

In general, the VEVs of the neutral Higgs fields can have non-trivial phases that induce the 
spontaneous CP violation. However, as shown in [13], there cannot be such spontaneous CP 
violation in the minimal PQ-NMSSM since the extrema with non-trivial phases are local 
maxima rather than minima. Accordingly, one can always set the VEVs of the neutral Higgs 
fields to be real. The equations of motion dVmggs/dH'^ = dVmggs/ dH^ = dVmggs/dS = 
at the vacuum are 



ml + /ie\ + + ^l±^{vl - vj) - 6eff/tan/3 = 0, (25) 



+ /ie'ff + X'vl + ^^^{v\ - vl) - Ks tan /3 = 0, (26) 



vs [ml + A^^^] +is- XAxVuVd = 0, (27) 
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where tan/3 = Vu/vd, f ^ = + f^, /Ucfr = At^s, and 6cfr = yUcffAA + Ayu|. These three 
minimization equations can also be cast into the following form, 

sin 2/3 = 3 ^^"^ „ r:-^ , (28) 

.2 



1 o mfr — mi, tan /3 „ 
2^^~ tan2/3-l 



o^l = ^^rrrrf^^^ - , (29) 



AAA^;2sin2/3-2t5 
= 2(m| + AV) ' 

where v ~ 174 GeV. In general, there can be false vacua that do not satisfy the proper 
EW vacuum conditions, i.e. v"^ ~ (174 GeV)^, v^, vs 7^ 0. For the EW vacuum to be 
stable, the false vacua should not be deeper than the EW vacuum or distant enough to take 
a longer time to decay than the age of the Universe. Finding such conditions, however, 
would be in need of systematic studies as carried out in [TSj, which consider the case of the 
Zs-invariant NMSSM. Although such extensive works in the case of the PQ-NMSSM are 
beyond the scope of this paper, we here leave a comment on the simplest condition for the 
false vacuum with Vu = = 0, vs 7^ o|3 For this false vacuum, the minimum value of the 
potential is given as 

^alse.min — 9" "I" l^S' (31) 

mi. 

On the other hand, the potential value at the EW minimum is 

,,,misin^2/3 micos^2/3 2ucfr Inn . 

^True,min = "A^ ' + + ^^^l/i^^ + /i^. (32) 

with g"^ = {gl + 5'|)/2. And, if one imposes the condition that Vrruc.min < Vpaisc.min, it can 
be shown that this is always satisfied for \ts\ ^ |AA|f^. For ts around the EW scale or 
positive value, l^Faise,min can be much deeper than Vrruc.min depending on the parameter 
values, or m| can be negative that it can break the stability of the vacuum. Therefore, in 
the following sections for the phenomenological analysis, we confine the parameter space 
to the large —ts ~ (9(TeV)^. The large negative ts value will make the singlet scalar 
heavy enough to be decoupled as can be seen in f l27|l . For the small ts < O(lOOGeV)^ 
region, a light singlet-like scalar or pseudoscalar in the Higgs sector still could give an 
interesting and different phenomenology [S]. However, more complicated investigation of 
the parameter space considering the above vacuum stability conditions should be performed 
for the region, which is beyond the scope of the present work. 

We now expand the Higgs scalar potential around the EW vacuum. By collecting 
quadratic terms and eliminating the soft mass terms through the minimization equations. 



""^For other false vacua, numerical analyses would be required even along simple field directions. 
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we get, for the CP-even Higgs fields in the {HdR, Hur, Sr) basis, 

(m| cos^ /3 + sin^ /3 (2A^f ^ — — m|) sin /3 cos /3 Xv{2ficscos P — Axsin (]y 
■m% sin^ P + cos^ /3 \v{2^cs sin P — Ax cos /3) 
m| + X^v'^ 

(33) 

where = 26efr/ sin 2/3. Upon diagonahzing the matrix, one eventually obtains the 
following tree-level lightest Higgs boson mass in the limit of ^ mz and ms ^ /iefr, A\. 

ml,^^ ^ ml cos' 2/3 + X\' (sin' 2/3 - (2^efr - sin2/3)^ \ ^ 

The second and third terms proportional to come from the NMSSM Yukawa coupling 
XSHuHfi in the superpotential, which can significantly enhance the tee-level lightest Higgs 
boson mass compared to the MSSM. The second term is due to the doublet scalar quartic 
coupling \H^H^\ and is sensitively becoming small for the large tan /3 (sin 2/3 ~ 2/ tan /3), 
while the third term is from the doublet-singlet mixing and A^|S'p|i/,°j;p terms in the scalar 
potential. 

For the CP-odd Higgs fields, 

(m\siv? (3 m^sin/3cos/3 At;y4Asin/3\ 
m^cos2/3 XvAxcosp\ (35) 
m| + X'v'/ 

in the {H^i, H^i, Sj) basis. After dropping the Goldstone mode obtained by rotating the 
upper 2x2 matrix, the following mass matrix appears in the {A, Si) basis. 

The corresponding mass eigenvalues are 

mX^Autree = ^(^^1 + ^1 + X\' T sJ (m| -m\ + A2i;2)2 + AX'v^Al^ . 

For the charged Higgs fields. 



(37) 



K ai , 1 1 si i\ ( cos^/3 cos/3sin/3\ 
M\ = {ml + - X'v') [ .20 ] 38 

\cos/3sm/3 sm /3 / 

in the (-ffj", H^^*) basis. By rotating the matrix by the angle 7r/2 — /3, one can obtain the 
Goldstone mode and a mass eigenstate with the mass eigenvalue, 

"^f/±,tree = m\ + m'^ - X%' . (39) 



All the above masses will receive important loop corrections (see Appendix C in [3]). 
For the analysis in Sees. |3] and O we include the loop contributions as implemented in the 
Nmssmtools 3.1.0 



2.3.2 Neutralino sector 



The singlet superfield 5* in the NMSSM also significantly changes the neutralino sector 
compared with the MSSM. The additional singlino field S mixes with the neutral Higgsinos 
H^, and the gauginos Ai, X^, producing a symmetric 5x5 mass matrix M^o, 



/Ml 





Ms 



-giVd/ V2 
g2Vd/ V2 




giVu/ V2 
-g2Vu/ V2 





V 










(40) 



in the basis (— iAi, —iX^, H^, -ff^, S). The diagonalization of this mass matrix and re- 
sulting mixing matrices are computed in Appendix |Al An important point to note is that 
there is no SUSY mass term for the singlino S in the minimal PQ-NMSSM. The singlino- 
like neutralino mass is induced only by mixing, and thus making the corresponding mass 
eigenvalue generically quite small. The lightest neutralino mass appears to be 



-2 ificsNisNu + Xv cos f3NuNi5 + A^; sin /JiVisiVis) 



+V2v (giNn 

x^v' r 

sin 2/3 



X^V^ 



-Ni3 COS (3 + Nu sin (3) + M^N^^ + M2N- 



12 



sin 2/3 



gjv^ 



2/ieffMi 2/ieffM2 



(41) 



when Mi,M2 ~ fica ^ Xv. Here, A'^ij's denote the neutralino mixing components given 
in Appendix |X1 The larger tan /3 makes the lightest neutralino mass smaller. Moreover, 
one can see that the lightest neutralino becomes lighter as fies increases. It can also be 
easily checked that the mass becomes zero as fics vanishes if we consider decoupling of 
the gauginos, i.e.. Mi, M2 ^ /ieg, Xv. This observation implies that there must exist a 
maximum value of the lightest neutralino mass at a certain value of fics for a fixed value of 
Xv in decoupling limit of the gauginos. We can find that the upper bound of the lightest 
neutralino mass is given by 



m^o = Xv cos /3 at = Xv sin (3 



(42) 



for tan (3 > 1, which gives m^o < 85 GeV for A = 0.7. 



^ The current version of the Nmssmtools implements only the Za-invariant NMSSM. We used modified 
codes for adapting to the case of the PQ-NMSSM. 
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3 Higgs phenomenology of the PQ-NMSSM 



In this section, we discuss the Higgs phenomenology of the PQ-NMSSM in the hght of 
the recent ATLAS and CMS discovery of a 125 GeV Higgs boson. We will show viable 
parameter spaces that can give the presumed Higgs boson mass without conflicting any 
existing phenomenological constraints, and discuss their feasibility from the point of view 
of naturalness. 

As discussed in Sec. 12.3. H the relevant parameters for the Higgs sector are A, tan/5, 
fJ'cS, nT'A (or Ax), fi's, and ts- For we take a weak scale value of 0(100 GeV)^ since it 
affects the Higgs sector only via the bes (= A^cfi^A + A/i|) term so that the variations of 
fies or niAiAx) include its effect. For and ts, we select their proper values that can 
give m| ^ m\ 3> m| in order to make the tree-level Higgs boson mass as large as possible 
according to fl34|l FI A large value of A is needed to increase the Higgs boson mass through 
the specific contribution of the NMSSM. However, an investigation into the scale of the 
perturbativity breaking by A driven by renormalization group running should be preceded 
to avoid the scale being lower than the PQ scale of 10^ - 10^^ GeV. It can be shown that 
the perturbativity breaking scale decreases rapidly as tan /3 becomes small. In Table [H 
we show the perturbativity breaking scales in the small tan /3 region for a fixed value of 
A = 0.7, which is a marginal choice for the PQ-NMSSM to be viable in the small tan (3 
region. 



Table 1: Perturbativity breaking scale A for tan/3 < 2. 



tan P 


1.1 


1.2 


1.3 


1.4 


1.5 


1.7 1.9 


A (GeV) 


lO^o 




10^2 


10^3 


10^4 





One of the most important features of the PQ-NMSSM is that the lightest neutralino 
is relatively lighter than that of the MSSM. Such a light neutralino can raise conflicts 
with several phenomenological constraints. Furthermore, since the neutralino sector shares 
some parameters with the Higgs sector, the constraints can also place serious restrictions 
on the viable parameter space for the 125 GeV Higgs signals. As for the neutralino sec- 
tor, the relevant parameters are A, tan /3, fi^s. Mi, and M2 as can be seen in Sec. 12.3.21 
Among them, the overlapping parameters with the Higgs sector are A, tan/3, and Hes- The 
phenomenological constraints on the light neutralino come from 

I. Z invisible decay {Z — )• XiXi); 

The region of Tn\ ~ will also be considered in the subsequent discussion with regard to the 
suppression of the ratio of the Higgs invisible decay. However, it will turn out to be unfavorable to the 
constraint on the charged Higgs boson. 
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II- XiX2 production at LEP II (e"^e — )■ XiX^); 
III. Higgs invisible decay {h — )■ XiXi)- 

These constraints impose additional restrictions on the parameter space of tan/3 and /ies- 
In particular, the constraints I and III motivate us to consider two scenarios according to 
the value of tan/3. For small tan/3 < 1.6, the lightest neutralino turns out to be heavier 
than the half of the 125 GeV Higgs mass so that the constraints I and III are satisfied by 
kinematics. For large tan /3 > 1.6, however, it should be considered a non-trivial mechanism 
of suppressing the Higgs invisible decay. 

For our study with specific choices of the usual MSSM parameters, the gluino mass 
parameter M3 is taken to be 1 TeV, which is above the region ruled out by current search 
results |T7], while being able to be accessed in the LHC experiment. For the same purpose, 
we set the first two generation squarks at ~ 2 TeV. In regard to the remaining gaugino 
mass parameters. Mi and M2, we consider two kinds of scenarios, the grand unification 
theory (GUT) relation of the gaugino masses and the unified gaugino masses in the TeV 
scale mirage mediation |18] . 

GUT : 6M1 = 3M2 = M3 = 1 TeV, (43) 
TeV mirage : Mi = M2 = M3 = 1 TeV. (44) 

In the case of the GUT relation, the lightest neutralino becomes generically lighter than 
in the case of decoupled heavy EW gauginos since the light EW gauginos can have a mix- 
ing with the Higgsino and the singlino. Furthermore, the mixing can deliver important 
phenomenological consequences through affecting the coupling of Higgs to the lightest neu- 
tralino pair as to be discussed in the following subsections. On the other hand, in the case 
of the TeV mirage, the EW gauginos are heavy enough so that the relation f H2]) holds to 
give heavier lightest neutralino, and the lightest neutralino is almost composed of Higgsino 
and singlino O 

In Subsecs. 13. H and |3.2[ we review the constraints I, II, and HI in the case of the light 
neutralino. We then proceed to analyze two scenarios depending on tan /3 in Subsecs. 13.3.11 
and [2321 

3.1 LEP II constraints on the light neutraUno 

In the PQ-NMSSM, not only the lightest singlino-like neutralino Xi but also the next- 
to-lightest neutralino X2 can be light enough to get constrained by the Z invisible decay 

''Here we concentrate on the minimal PQ-NMSSM. Phenomenological implication of the TeV-scale 
mirage mediation in the singlet extension of the MSSM (|17p was discussed in [T^ . 
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and/or the direct production in the e^e scattering of the LEP II experiment through the 
processes shown in Fig. [H 

The Z boson decay rate to a lightest-neutralino pair is given by 



9l (N' 
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(45) 



where A'^ij's are neutralino mixing components shown explicitly in Appendix |X1 From the 
constraint on the Z invisible decay, Finv < 3 MeV [20] , we find that 

(46) 



1^13 - Nil < 0.13. 



In the LEP II experiment, the next-to-lightest neutralino can be produced in associ- 
ation with the lightest neutralino as in Fig. [1] If one assumes that the sleptons are much 
heavier than the Z boson, the slepton exchange diagrams can be neglected in the LEP II 
experiment. Then, the production cross section is given as 



1 
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The production cross section is given upper bounds depending on the masses of the neu- 
tralinos by the OPAL analysis result [21] . In the case that the lightest neutralino is nearly 
massless, the most conservative bound is given as follows. 



a(e+e-^X?X2)<10fb 



(48) 
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with Br(x2 — )■ Zxi) = 1- However, such an upper hmit of the production cross section can 
be moderated in the case that the hghtest neutrahno is heavier. This will be discussed in 
more detail in Subsec. I3.3.1I 



3.2 Higgs invisible decay 

The light neutralino states can also induce the invisible decay of the Higgs boson. If there 
exists substantial invisible decay ratio of the Higgs boson, the visible decay modes such 
/i — 77 or /i — WW/ZZ will be significantly reduced. As this would make it difficult to 
explain the recent discovery of the 125 GeV Higgs boson signal in the ATLAS and CMS 
experiments, it is of particular importance to consider the Higgs invisible decay in the 
plausible parameter space. 

The invisible decay of the lightest CP-even Higgs boson is determined by the Higgs- 
neutralino coupling shown in Appendix [Cl 

g,^^ = A(^^^n5 + ^i2Hjf + ^i3Hj?) + |(^nHlf-5i2H^)-|(5nH^-^i2Hg^), (49) 

where U^j = NiaNjb + NibNja- Here, Sab and Nij are the CP-even Higgs mixing matrix and 
the neutralino mixing matrix, respectively. As analyzed in Appendix |A] including the first, 
second, and dominant third-order contributions, the mixing components in the neutralino 
mass matrix are 



QiXv"^ cos 2 13 
•11 — 

g2Xv'^cos2(3 
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Xv COS 13 Xv^ f gl \ np ■ p 3A^t;^ 



(50) 
(51) 



Xv sin /3 



cos 2/3 sin /3 ^ ^ sin 2/3 sin /3, (52) 

cff 

(53) 



iVi5 ^1. (54) 

We here include only the leading term for each component except Ni^, for which the second 
and third terms in ( l52l) can be sizable when tan/3 > [n/XvY. By putting these terms into 
POj) . one can find ghx\x° coupling between the lightest CP-even Higgs and a lightest- 

neutralino pair. In the limiting case of vanishing doublet-singlet mixing in the Higgs sector, 
i.e., S'13 = 0, the coupling is given by 



V2X'^v ( ^m%r 3m2 



/ieff V 4M2Veff M2/ieff 



tan B , , 

(55) 
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where c = 511/512 and g2 — 2gi. It was assumed that the GUT relation for the gaugino 
masses, M2 = 2Mi, and apphed a crude approximation, sin /3 ^ 1, cos/3 ~ l/tan/3, 
sin 2/3 ~ 2/ tan/3 and cos 2/3 ^ 1 for tan /^ ^ 3 to derive the above relation. Note that the 
coefficient c corresponds to — tana (— 7r/2 < a < 0) in the MSSM limit, where a is the 
CP-even Higgs mixing angle in the MSSM. For making this coupling vanishing, we find a 
relation, 
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(56) 



2M2/ieff 

For M2/icff ^ "^vv' t'^^ coefficients of c and are always positive, so the larger c requires 
the larger value of tan /3. The effect of the Higgs mixing c will be discussed in Subsec. 13.3.21 



3.3 Phenomenology depending on tan f3 
3.3.1 1 < tan /3 < 2: heavy neutralino scenario 

We now discuss the possibility for the small 1 < tan /3 < 2 region. Such small tan /3 makes 
singlino-Higgsino mixing large and thus the lightest neutralino heavy. One can obtain 
the mass of the lightest neutralino larger than mh/2 in order to forbid kinematically the 
invisible decays of the Z boson and the Higgs boson. In this case, however, it may concern 
the x5x2 P^i'^ production as the next-to-lightest neutralino can also be light enough. As 
discussed in Sec. 13.11 the most conservative bound for such a neutralino-pair production 
cross section was set by the OPAL result. The upper bound of the cross section varies 
according to mass values of the neutralinos. It becomes as much as 70 fb in the parameter 
space of our interest, m^o > mh/2 ~ 63 GeV and 120 < m^o < 140 GeV. 

The numerical results for the GUT and TeV mirage relations of gaugino masses are 
shown in Figs. [2] and [3l respectively. In the case of the GUT relation, light EW gauginos 
are mixed with the singlino, and the lightest neutralino becomes lighter. Moreover, the 
next-to-lightest neutralino also becomes lighter due to the Higgsino-gaugino mixing, which 
increases the X1X2 production in the LEP II experiment. Hence, only small region of the 
parameter space, tan/3 < 1.2 and 120 GeV < /ics ^ 130 GeV, is allowed. On the other 
hand, in the case of the TeV mirage, the EW gauginos are heavy so that the lightest 
and the next-to-lightest neutralinos are mostly composed of the singlino and the Higgsino. 
Consequently, the masses of the lightest two neutralinos are heavier than those of the 
GUT relation case, and thereby avoiding the LEP II constraint for the X1X2 production for 
broader parameter region, tan/3 < 1.5 and 130 GeV < fics ^ 170 GeV. 

This small tan /3 scenario has excellent features in the naturalness point of view. The 
large loop corrections by the stops are not necessary to raise the Higgs boson mass up to 
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Figure 2: Plot for small tan/3. Here we set 6Mi = 3M2 = Afa = 1 TeV, tua = 350 GeV, 
"^Qs ~ ~ GeV, At = 0, and ^5 = —7 x 10'' GeV^. Black curves denote Higgs mass in 
GeV, blue curves denote m~o. in GeV, and green dashed curves denote tu^o. The gray-shaded 
region is excluded by the OPAL [21]. 

125 GeV as the tree-level Higgs boson mass can be raised enough by the A-proportional 
contribution by virtue of the small tan /3. This ameliorates the fine-tuning problem from the 
stop sector. Moreover, all massive soft parameters related to the EW symmetry breaking 
are of order of 100 GeV. When tan/3 = 1.3 and fies = 135 GeV in Fig. |3l for example, the 
masses related to the EW symmetry breaking are given by 

rriH^ ~ rriH^ ~ 200 GeV, ms ~ 600 GeV, mj ~ 500 GeV. (57) 

Therefore, the EW symmetry breaking condition can be satisfied up to 5% parameter 
tuning. 

Concerning the cosmology in this scenario, if the LSP is the lightest neutralino, the 
substantial Higgsino mixing in the LSP becomes the cause of the similar cosmological 
features as the Higgsino dark matter of the MSSM. Still, due to the large Higgsino-singlino- 
Higgs coupling from the XSHuHd term in the superpotential, the direct detection cross 
section can be rather large. On the other hand, it is expected that the amount of missing 
energy in the LHC experiment will not be much different from that of the MSSM since 
the rriy^o is slightly larger than mh/2. The existence of the light third-generation squarks 
will lead to the top or bottom-rich signal. The dark matter and the collider signatures 
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Figure 3: Plot for small tan (3. Here we set Mi = M2 = M3 = 1 TeV, tha = 350 GeV, = 
m-^c = 500 GeV, = 0, and S^s = — 7 x 10^ GeV'^. Black curves denote Higgs mass in GeV, 
blue curves denote m^o. in GeV, and green dashed curves denote m^o. The gray-shaded region 
is excluded by the OPAL [S]. 

at the LHC for the small tan /3 region will be discussed in more detail in Sees. H] and |5l 
respectively. 

3.3.2 Large tan /3: very light neutralino scenario 

For large values of tan/3 > 1.6, the lightest neutralino mass becomes smaller than mh/2, 
and the Higgs invisible decay mode to a lightest neutralino pair is open. In this region, 
it is necessary to examine conditions for suppressing the Higgs invisible decay. As argued 
in (!56|) . it is plausible to put the Higgs mixing parameter c as small as possible in order 
to make the Higgs invisible decay rate vanishing, while keeping relatively small tan /3 by 
which the quartic coupling A can substantially raise the tree-level Higgs boson mass. As 
an extreme example, we consider vanishing Higgs mixing parameter c. To achieve this for 
relatively small tan/3, we set Ai'su = implying 



(2A^t'^ — m\ — m|) cos /3 sin /3 = 0, 



(58) 



that is, 



2 o \ 2 2 2 

= 2A V — rriz 
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(59) 
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In this case, the hghtest Higgs boson is mostly the up-type, Hu, for tan f3 > 2. Here, m\ 
is not the physical CP-odd Higgs boson mass, but a model parameter defined as m\ = 
26cff/ sin 2/3, which indeed corresponds to the CP-odd Higgs boson mass in the MSSM limit. 
Now, by putting c = in fl56|) . one can find 



and thus tan /3 3 for /^eff = ^2 = 300 GeV. Although one can suppress the Higgs invisible 
decay with tan /3 ~ 3 in this way, the small m\ in ( 159|1 can lead to the light charged Higgs 
boson. At tree level, the charged Higgs boson mass is given by 



leading to m/^± ^113 GeV for the parameter choice ( l59ll . Such a light charged Higgs boson 
is excluded by the recent ATLAS search in the decays of top quarks |22]. More generally, 
the charged Higgs boson mass smaller than 150 GeV and tan/3 < 4 has been excluded 
under the assumption of Bt{H~^ — t~^u) = 1. Fig. H] shows the result of our calculation of 
the Higgs invisible decay ratio depending on ruA and tan/3 and the current LHC limit. It 
should be mentioned that in the horizontal axis is the input parameter at stop mass 
scale, so it is different from in (!58|) that is the value at the weak scale. In addition, since 
the result includes loop corrections, it is slightly different from what is expected from the 
tree-level estimation. In Fig. HI one can find the region of tan /^ > 4 and rriA ~ 130 GeV in 
which the Higgs invisible decay ratio becomes vanishingly small. In this region, however, 
the production cross section and the decay branching fraction of the lightest Higgs boson 
become different from the SM-like one because of the CP-even Higgs mixing from the mass 
matrix fl33|) . In order to assess the viability of such non-SM-like Higgs scenario, we would 
need more complicated study for whole SUSY parameters and collider signatures. We will 
leave this for future work and focus on the SM-like Higgs phenomenology from now on. 

In the case of larger m^, the invisible decay branching fraction of the Higgs boson 
drastically increases so that larger tan /5 is needed for the Higgs invisible decay to vanish 
as was showed in (!56|) . For rriA ^ 200 GeV, the invisible branching fraction becomes nearly 
independent of rriA since such region corresponds to the MSSM decoupling limit, in which 
the lightest CP-even Higgs boson becomes SM-like, so we find c ^ — tan a ^ cot /3. In other 
words, the Higgs coupling to a neutralino pair f l55|) depends only on tan p. In this region, 
we need tan /3 > 9 for BT{h — )■ XiXi) < 0.1. It should be noted that this result is obtained 
with the GUT relation of gaugino masses, and if one considers the TeV mirage relation, the 
corresponding tan/3 value becomes larger. For such large tan/3, the NMSSM feature of the 
sizable A contribution to the tree-level Higgs boson mass is lost since A^f ^ sin^ 2/3 in flM|) 
is already small compared to the Z boson contribution m| cos^ 2/3. Thus, it is required to 
have large loop corrections from the stop sector as in the MSSM. 
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Figure 4: Contours of the branching fraction of the Higgs invisible decay for ^Ucfr and the GUT 
relation of gaugino masses. In the left side of red dashed region, charged Higgs mass is smaller 
than 150 GeV so that tan/3 ^ 4 is excluded by the charged Higgs search experiment. 

For the comparison between cases of the GUT and TeV mirage scenarios, we show the 
Higgs couplings to a neutralino pair and corresponding branching fractions for both cases in 
Figs. [5] and |6l In the case of the GUT relation, the EW gauginos are light (Mi = 160 GeV 
and M2 = 330 GeV), so the coupling vanishes when tan /3 ~ 12 as can be seen in Fig. [Sl(a). 
This is consistent with the relation fISB]) . Therefore, the branching fraction of the Higgs 
invisible decay is very small near tan/3 ~ 12 as shown in Fig. [5](b). On the other hand, in 
the case of TeV mirage relation, the gauginos are relatively heavy (Mi = M2 = 1 TeV), 
thus (156|) can be satisfied only for very large values of tan (3. In Fig. [6](a), we cannot see the 
solution of (15^ in the range of 1 < tan (3 < 55. Instead, the size of the Higgs coupling to 
a neutralino pair becomes smaller as tan /3 becomes larger. Hence, nearly vanishing Higgs 
invisible decay can be attained when tan /3 > 20 as shown in Fig. [6](b). 

So far, we have seen that rather large tan/3 (9 < tan/3 < 18 for the GUT and 
tan/3 > 20 for the TeV mirage relation) is needed to reduce the Higgs invisible decay 
branching fraction below 10% level. Even if we allow the branching fraction up to 50% 
level, one needs tan /3 ^ 6 for the GUT and tan /3 > 10 for the TeV mirage relation. 
This suppresses the NMSSM contribution to the tree-level Higgs boson mass and there 
is no difference from the MSSM in the naturalness point of view. However, an important 
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Figure 5: (a) The lightest Higgs couphng to the hghtest neutralino pair, (b) Branching fraction 
of the Higgs invisible decay. In these plots, we use tua = 500 GeV, fieS = 400 GeV and the GUT 
relation of gaugino masses. 
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Figure 6: (a) The lightest Higgs coupling to the lightest neutralino pair, (b) Branching fraction 
of the Higgs invisible decay. In these plots, we use uia = 500 GeV, fics = 400 GeV and the TeV 
mirage relation of gaugino masses. 



difference arises due to the existence of the very light neutralino. As one can see in fj^Tj) . the 
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lightest neutralino mass is generally smaller than 15 GeV for tan /3 > 5 leading to non-trivial 
implications in cosmology and collider signatures. If such a light singlino-like neutralino is 
the LSP, its annihilation cross section is very small, and thereby overproducing the dark 
matter density. To avoid this problem, we need a non-standard cosmological history, which 
will be discussed in Sec. HI On the other hand, the singlino-like LSP modifies the decay 
topology of the supersymmetric particles. This will be discussed in Sec. |5l 

4 Dark Matter 

Let us now discuss the dark matter cosmology of the singlino-like LSP for two regions of 
small and large tan /3 that survived the various constraints from particle phenomenology. 
For each region, the standard LSP dark matter obtained by thermal freeze-out has a dif- 
ficulty in satisfying the cosmological density and/or the direct detection bound. Instead, 
the non-standard cosmology driven by the axion supermultiplet, the PQ sector, can pro- 
vide a viable range of the parameter space. In the following discussions, we assume that 
R-parity is conserved and the lightest neutralino, mostly singlino, is the LSP disregarding 
a possibility of lighter axino or gravitino being the LSP. 

4.1 Small tan/3 region (1 < tan/3 < 2) 

In the small tan /3 region, the LSP mass m^o is around 50 — 70 GeV. As studied in the 
context of the nMSSM [23], both s-channel Z boson exchange {xiXi Z ^ ff) and s- 
channel Higgs boson exchange {xiXi h ff) are equally important for the neutralino 
annihilation cross- sect ion. The effect of the Higgs exchange is more important than in 
the case of the MSSM since the h-XrXi coupling (!66|) is enhanced compared to that of the 
MSSM. The resulting relic abundance fl-^h"^ is represented in the left panels of Figs.[7]and[8] 
for the GUT scale unified gaugino masses and the TeV scale mirage mediation, respectively. 
The superscript "TH" implies the dark matter density produced from thermal freeze-out. 
In each figure, the narrow contour lines of around 10~^ — 10~^ for smaller tan/3 

correspond to the Higgs resonance region, m^o ^ rah/2. As tan /3 increases on the right- 
hand side of the Higgs resonance, m^o < m/i/2, m^o approaches mzjl and thus fi™/;,^ 
decreases again due to the Z resonance. As tan /3 decreases on the left-hand side of the 
Higgs resonance, m^o > mhll^ the neutralino relic abundance increases and can reach the 
correct amount of the dark matter, ViJ^h? ^ 0.11. 

Right panels of Figs. [7] and [8] represent the effective spin-independent elastic scattering 
cross-section between the nucleon and the dark matter in units of 10~^ pb. For the direct 
detection bound on the scattering cross-section [S] , the standard local dark matter density 
is assumed, that is, nx = Pbm/^xi — ^■'^ cm~^(GeV/m^o) corresponding to ^dm^^ ~ 0.11. 
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Figure 7: Plot for low tan /3 with 6Mi = 3M2 = M3 = 1 TeV. Other parameters are the same as 
those of Fig. [2j (a) Black curves denote the thermal relic density of the neutralino LSP, r2™/i^. 
(b) Red curves denote the central values of the nucleonic scattering cross-section of the lightest 
neutralino, (i7™/riDM)<7gi in the unit of 10~^ pb. The magenta curve denotes 125 GeV Higgs 
mass, the purple curve denotes m~o = m/j/2, and the gray-shaded region is excluded by OPAL. 
For 50 GeV < m~o < 70 GeV, the XENONIOO rules out the region (I^^o^/J^dm) cr^i^^o > 2 x 10"^ 
pb at 90% C.L. p]. 



However, in our parameter region where the thermal relic abundance of the neutralino dark 

matter is smaller than QDuh"^ ~ 0.11, the neutralino number density around the Earth 

deceases by the factor of f2™/f2DM- Thus, we show the effective cross-section between the 

Xi 

nucleon and the assumed dark matter, {Q-^/Qdm) x (^si\°^ where the spin-independent 
elastic scattering cross-section between the nucleon and the neutralino crsi,x1 given by 

^si.xO = - ] fl (62) 

Here, represents the nucleon [p, n) and the constant is 

^ = E C (-) + ^ E C (^) . (63) 

q=u,d,s ^ ^ q=c,b,t ^ ' 



where rriq denotes the quark mass, is the coupling constant of the quark level effective 
Lagrangian term £eff = ctgXiXi^Q'j and the hadronic matrix elements are given by m^f'^^ = 
{N\mqqq\N) and f!j^^ = 1 — ^^^^ ^ ^ [21] • The numerical values of the matrix element 
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Figure 8: Plot for low tan /3 with Mi = M2 = M3 = 1 TeV. Other parameters are the same 
as those of Fig. [2l (a) Black cm'ves denote the relic amount of the neutralino LSP, r2™/i^. 
(b) Red curves denote the central values of the nucleonic scattering cross-section of the lightest 
neutralino, {Q,™/Qdm) <7gj in the unit of 10~^ pb. The magenta curve denotes 125 GeV Higgs 
mass, the purple curve denotes m~o = m/j/2, and the gray-shaded region is excluded by OPAL. 
For 50 GeV < m~o < 70 GeV, the XENONIOO rules out the region (J7To"/J1dm)o-si,^o > 2 x 10"^ 
pb at 90% C.L. §]. 



are determined in [25] 

/^P) = 0.020 ± 0.004, = 0.026 ± 0.005, /^^j = 0.118 ± 0.062, 

/^^) = 0.014 ± 0.003, = 0.036 ± 0.008, = 0.118 ± 0.062. (64) 

In our case, the coupling constant aq is dominantly given by t-channel Higgs exchange 
diagram, 

ml \mwsmpj 

for the up-type quarks. For the down-type quark, aq is obtained by the appropriate re- 
placements {Sii —J- 5*21, sin/3 — )■ cos/3). The Higgs-LSP-LSP coupling Qhxlxl given by 

9hx1x\ = -92{Ni2 - tan6'H/iViij {SuNn - S'laA^u) 

+ V^A(Si3iVl3iVi4 + iVl5(^12iVl3 + ^lliVl4)). (66) 



The first line in (!66|) is the same as in the MSSM, and the second line comes from the 
superpotential term, XSHuHd- Due to the sizable Higgsino component of the LSP, the 
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second term dominates (jf^^o^o. The recent XENON 100 data puts an upper limit on the 
effective spin-independent LSP-nucleon cross-section around 2 x 10~^ pb for the mass range 
of 50-70 GeV at 90% C.L. [9]. Note that the region that gives nl^h^ = 0.11 is far above the 
XENONIOO bound. However, for the points on the 125 GeV Higgs hne close to the Higgs 
resonance point, the XENONIOO bound can be avoided because the neutrahno thermal 
rehc density can be much below the observed dark matter density. More specifically, we 
find that the mass difference between m^o and mh/2 should be smaller than (9(0.1) GeV, 

r«x?-^^0.1GeV. (67) 

Recall that, since the neutralino contribution to the dark matter density is small, the major 
component of the dark matter density can come from the axion for vpQ ~ 10^^ GeV in our 
scenario, 

"»"'"°'2Kio5¥liv)"°<«'>- 

where 6 is the initial misalignment angle, typically of order one. 

Let us note that the PQ-NMSSM has a late-time decaying saxion field s, a CP-even 
partner of the axion a, which can dilute away the LSP relic density calculated previously. 



If the PQ sector is stabilized by SUSY breaking effects as discussed in Sec. 12.11 the saxion 
typically gets a mass of the order of the weak scale, and its interaction with other particles 
are suppressed by l/upq, so it has a long lifetime. It is natural to have a period in the 
early universe during which the energy density of the universe is dominated by the coherent 
oscillation of s or vacuum energy of the PQ sector field. At the end of this period (t ~ 1 /Tg 
where is the total decay rate of s), the saxion will decay to produce radiation reheating 
the Universe at the temperature Trh obtained from the relation Ps(^rh) — Pr(rRH)- The 
reheat temperature Trh is given by 

. 200 MeV (-p^) (-) (_) (_— ), ,69) 

where Br^ is the branching fraction of the saxion decay into the axion pair having the rate 
rs_s.aa = ''^s/(647rUpQ) = PsBra [26]. For vpQ = 10^^ GeV, the reheat temperature is much 



smaller than the freeze-out temperature of the LSP (Trh ^ ^/ — '^x?/^^ ~ GeV), and 
thus the depleted thermal LSP population cannot be regenerated. Furthermore, the non- 
thermal production of the LSP from the saxion decay can also be forbidden if the saxion 
mass rris is taken to be smaller than 2m^o. Therefore, the stringent constraint (167|) on the 
LSP mass can be relaxed. 



4.2 Large tan f3 region (5 < tan f3) 

In the large tan /3 region, we found that the lightest neutralino mass is below 0{5 GeV) from 
the consideration of the Higgs invisible decay. Since m^o is far from the Higgs resonance 
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region, the LSPs are annihilated into the SM fermion pairs dominantly by the s-channel Z 
boson exchange, Z ^ f f. In the mass range of 1 GeV < m^o < 5 GeV, the LSP 

rehc abundance is given by 

Q^fh' = 1.07 X 10^ (^7^^) (^"^^ dTg,{TY/'{a,^,,v)T GeV'^ (70) 

10 \ / 4m|o V /2 GeV\ ^ / 0.04 



.S.(Tf)J \ ml j \ mj,, j \\Nl-NU 

where Tf is the freeze-out temperature of Xi? = m^o/Tf is around 8 ~ 10 (15) for 
m^o ~ 1 GeV (5 GeV), and g^Tf) ~ 10 for T/ ~ 100 MeV. The order-one constant 
7jjjro basically counts the number of fermions lighter than the LSP. The neutralino mixing 
elements f lS^ and f lS^ give 

,..2 ..2, AVsin2/3 / m^o N /tan/3\ /500 GeV\ 

|/V?,-A'f,N^^.0.04(^)(^)(^_) (71) 

for large tan/3. Recall that liVfg — N^^\ is bounded by 0.13 as in ( H6l) . Therefore, in order 
to avoid the overclosure dark matter density, the neutralino abundance has to be depleted 
by the factor of A ~ 10^ - 10^ for m^o ~ 1 - 5 GeV. 

In the large tan/3 region, m^o is sensitive to A while the Higgs mass is not, and can 
be much smaller than 1 GeV if we allow for A smaller than the nominal choice 0.7. For 
instance, one finds m^o ~ 10 MeV for A ~ 0.1. In this case, the freeze out temperature Tf 
is larger than m^o and becomes a few GeV. The relic abundance is then given by 

"?'"^--io^(3^)(-|L). 

Thus, even larger dilution factor A ~ 10^ is needed for m^o = 0(10 MeV). 

It is amusing to note that the dilution mechanism by the saxion field discussed in 
the previous subsection can successfully deplete the LSP thermal abundance as well as 
produce the right amount of the non-thermal LSP relic density. Since the singlino LSP is 
light enough to be produced by the saxion decay, we need to check how sizable amount of 
the LSP dark matter can be produced in this process. Such a non-thermal production is 
controlled by the ratio m^o/rris as the s-XrXi coupling is proportional to m^o/vpQ whereas 
the s-a-a coupling is proportional to ms/fpg. We here consider the case ~ 100 GeV, 
then the dominant saxion decay mode can be s — )■ fefe whose coupling is proportional 
to Csz/mfe tan/3/fpQ, where Csh is an order-one parameter representing the saxion-Higgs 
mixing. For tan /3 > 25, the s-b-b coupling can be larger than the s-a-a coupling allowing 
for Bra < 0.1. Now, one can find the relic abundance of non-thermally produced LSP 
through the decay s — XiXi as follows. 

aS.™/.= = 0.1-,,,.,-. C-^) (^) , (73) 

^° ^'^1^1 vsoMeV/ V yVioMevy' ^ ^ 
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where 7s^5x? is the order-one constant controlhng the s-XrXi coupling. The above relation 
shows that the light singlino LSP abundance can be in the right range for m^o ~ 30 MeV 
(obtainable for, e.g., A ~ 0.07, tan/3 = 25, fi^s = 300 GeV), m, ~ 100 GeV, and vpq ~ lO^^ 
GeV. Since we used the PQ symmetry breaking scale vpq larger than 10^^ GeV, we also 
have to consider the axion dark matter abundance. The entropy dumping from the saxion 
decay at Trh = O{10) MeV can dilute the axion relic density as well. In this case, the 
axion dark matter at low decay temperature is given by [26] 

(I 0-98 1 Q 

f^(^^)'] ( )"f™^fevV (^^ (74) 

V 10 VlOMeVy J VlOi^GeV/ V ^QCD / ^ ^ ^ 

which is negligible. 

5 Collider signature at the LHC 

In this section, we study the collider signature of the PQ-NMSSM taking some benchmark 
parameter points in the case of the TeV scale mirage relation of gaugino masses discussed 
in the previous sections. The mass spectra and branching ratios of the sparticles and Higgs 
bosons have been calculated with the modified codes of Nmssmtools, and are given in 
Tables [2] and [3] in the small tan /3 scenario, and in Tables [Hand [5] in the large tan /3 scenario, 
respectively. Notice that the stop mass is chosen to be as light as ~ 500 GeV consistently 
with the 125 GeV Higgs boson for the small tan/? benchmark point, whereas it has to be 
heavy (~ 1 TeV) with a large mixing for the large tan/3 benchmark point. 

In the small tan/3 scenario, the main decay topologies are similar to those of the 
typical MSSM scenario with light stops and sbottoms due to the large mixing between the 
singlino and the Higgsino. On the other hand, in the large tan /3 scenario, the existence of 
a very light singlino LSP prevents the direct decay of the sparticles into the LSP and thus 
there appear additional decay steps compared to the similar MSSM decay processes. The 
main decay processes of the large tan /3 benchmark point are 

.g^t^t^ xfbt ^ XiWbt, 

• g^ht^ Xltt ^ x\Z (Hi) tt, 

• g ^hb ^ xftb ^ XiWtb, 

• g^bib^ xlbb XiZ (Hi) bb. 

Hereafter the SM-like Higgs boson h is denoted by Hi representing the lightest one among 
three CP-even Higgs bosons. As one can see from the above decay chains, the final states 
consist dominantly of multi-jet coming from the multi-top/bottom, and Hi/Z/W + missing 
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Table 2: Sparticle and Higgs masses (in GeV) for the small tan (3 benchmark point. 
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Table 3: Main decay modes for the sparticles and the lightest Higgs boson and their branching 
fractions in % for the small tan (3 benchmark point. 
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bi 
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99.1 







energy that may escape the early LHC searches on the channels with (two to four) jets 
and large missing energy [27] |^ Furthermore, the heavier neutralinos will eventually decay 
into the light singlino LSP, which makes the missing energy generically small as claimed 
recently in [30]. This will reduce the efficiency of the searches on the channels with large jet 
multiplicities and missing energy [3T||32] . Another notable feature is that the relatively light 
stops and sbottoms lead to top/bottom-rich signal events. The recent constraints on the 
stop and sbottom masses were set by the searches on the gluino-mediated production [33] 
and the direct production of stops and sbottoms [SI EH]. In our study, the gluino and 

^ This was also claimed in [551[ini in the context of Zs-invariant NMSSM. Here, we stress that the 
relatively light singlino-likc LSP can be naturally attained in the PQ-NMSSM. 



27 



Table 4: Sparticle and Higgs masses (in GeV) for the large tan/3 benchmark point. 
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Table 5: Main decay modes for the sparticles and the lightest Higgs boson and their branching 
fractions in % for the large tan /? benchmark point. 
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squarks are chosen be heavy enough to evade such search results. 

In order to study the collider phenomenology, we have generated Monte Carlo (MC) 
event samples for a proton-proton collision at 14 TeV by HerwigH — h 2.5.2 [H]. The 
generated event samples have been scaled to 10 fb~^ integrated luminosity. For the sake 
of the simple analysis, we consider only the leading-order cross sections calculated with 
HerwigH — h and the CTEQ6L [37] parton distribution functions (PDF). The generator- 
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level event samples have been further processed with the fast detector simulation program 
Delphes 2.0 [3H] using the ATLAS detector card. Jets are reconstructed using the anti- /ct 
jet clustering algorithm [39] with radius parameter of 0.4. Isolated electrons (muons) are 
required to have the transverse momentum pt > 20 (10) GeV and the pseudo-rapidity 
|?7| < 2.47 (2.4). In the recent version of Delphes, the missing transverse momenta fix 
are defined as the negative vector sum of the transverse momenta of all the calorimetric 
cells and muon candidates. On top of that, we resolve overlaps between jets with |?7| < 2.8 
and leptons by following the recent ATLAS analysis [31]. Jets lying within a distance 
AR = (Arjy + (A0)2 < 0.2 from an electron are discarded. Then, any lepton remaining 
within a distance AR < 0.4 from such a jet is discarded. Finally, all jets with |?7| > 2.8 are 
discarded. From now on, we will use only the remaining electrons, muons, and jets for the 
analysis. 

The dominant SM backgrounds for the signal with large jet multiplicities come from 
the top-pair process and W oi Z bosons produced in association with jets. The di-boson 
and the single-top processes can contribute to the backgrounds as well, but they are neg- 
ligible. For a simple study of the background effects, we here consider only the top-pair 
process. Since our study mainly concerns the signal events with a large number of jets, 
we generate fully hadronic, semi-leptonic, and fully leptonic ti events with up to two addi- 
tional partons in the matrix element using Alpgen 2.14 [3D] and CTEQ6L PDF sets. The 
parton showering to generate additional jets, and the fragmentation and hadronization are 
performed by HerwigH — \-. The MC samples have been processed with the Delphes to 
reconstruct jets and isolated leptons and adjust the detector effects. 

To suppress the backgrounds, we impose the basic event selection cuts as follows. 



(i) At least 6 jets with > 80 GeV, 

(ii) no isolated electron or muon, 

(iii) missing transverse energy > 80 GeV, 

(iv) St > 0.2. 



Here, is the transverse sphericity defined as 

5x - (75) 

where Ai and A2 are the eigenvalues of the 2x2 sphericity tensor = 'YlikPkiP^'' of the 
reconstructed objects [H]. This variable is known to be useful to suppress QCD events in 
which back-to-back configurations (iSx ~ 0) are dominated. Although such backgrounds are 
not considered here, this variable will be included in the analysis. In the ATLAS MC study 
for the inclusive SUSY search [H], the azimuthal angular separation cut, A0(jet — fix), is 
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/s = 14 TeV 



/s = 14 TeV 




Figure 9: The distributions of (a) the jet multiplicities and (b) the ^t/\/^t for the signal events 
and tt+ jets. 

imposed further to reduce the jet mis-measurement effect. However, it has been recently 
noted that this cut variable may lose the efficiency if the next-to-lightest neutralino X2 
boosted and the lightest Higgs boson Hi or Z boson decayed into bb aligned with the fix 
vector [29]. We also note that the cut on the missing transverse energy is required to reduce 
the fully hadronic ti background process, even though it should be relatively milder than 
the typical SUSY searches because of the small missing energy. 

We show the jet multiplicity distributions for the signals and backgrounds after apply- 
ing the basic selection cuts in Fig. [9](a). One can see that there are relatively more number 
of energetic jets in the large tan /3 scenario. This is because the additional decay step 
caused by the existence of the singlino LSP as well as the fact that the sparticle production 
modes are mostly consisted by a gluino-pair production as the squarks are too heavy to be 
produced. On the other hand, the tt and bb production cross sections are somehow sizable 
in the small tan /3 scenario. The relatively short decay steps result in the smaller number of 
jets. Fig. Mjo) shows the distribution of the ratio -^t/v^^Tj where i^x is the scalar sum of 
the transverse momenta of all jets with px > 50 GeV and |?7| < 2.8. This variable has been 
known to be useful to increase the performance of the missing energy reconstruction [32] . 
We require the condition of / V Ht being larger than 4 GeV^^^ in addition to the basic 
event selection cuts. By employing this cut variable on top of the basic selections cuts, 
the jet multiplicity distributions are shown in Fig. [TOT a). One can see that the background 
has been reduced, while the signals remain almost untouched in the region of the large jet 
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Figure 10: The distributions of (a) the jet multiphcities and (b) the M^s for the signal events and 
tt+ jets with imposing the / VHt cut. 

multiplicity. In Fig. [TOT b). we also show the distribution of the effective mass defined as 

Meff = ^T + 5^PT, (76) 

jets 

where the summation is over all jets with > 50 GeV and |?7| < 2.8 in the event. The 
Meff has been known to be a good variable for discriminating SUSY signal events from 
SM backgrounds |12]. In the literature, it was noted that the peak position of the M^s 
distribution has a strong correlation with the SUSY mass scale. In the case of the small 
tan/3 point, the Mgg distribution is peaked in the lower position than in the case of the 
large tan /3 point. This is because of a sizable amount of direct production rates of the light 
stops and sbottoms, whereas the signal events for the large tan /3 point come practically 
from the gluino-pair process. 

For a crude estimation of the signal significance, we show how the cross sections of the 
signals and the backgrounds change under each event selection cut in Table El Although 
we did not perform a more accurate estimation which requires higher-order cross sections, 
optimization of the cut values by the multivariate analysis techniques, and understanding 
the systematic uncertainties, we can expect the 5a discovery for both PQ-NMSSM scenarios 
with several fb~^ of the integrated luminosity. 
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Table 6: Cut flows of the signals and backgrounds in fb. 
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PQ-NMSSM 
small tan /3 large tan /3 


tt 


tt + 1 jet 


tt + 2 jets 


> 6 jets with > 80 GeV 


105.7 


57.3 


450.2 


1650.6 


2055.8 


Lepton veto 


51.8 


20.5 


359.0 


1259.8 


1528.9 


> 80 GeV 


44.7 


19.4 


29.7 


148.1 


219.1 


St > 0.2 


38.0 


16.6 


23.9 


119.6 


173.0 


^t/Z^t > 4 GeV^/2 


28.6 


13.8 


12.7 


58.8 


82.8 


Meff > 1000 GeV 


26.2 


13.8 


8.4 


45.9 


67.4 


Meff > 1500 GeV 


16.6 


12.4 


1.0 


11.5 


19.5 



6 Conclusions 

Motivated by the axion solution of the strong CP problem and the recent discovery of a 
125 GeV Higgs boson at the LHC, we investigated phenomenological consequences of the 
PQ symmetry realized in the context of NMSSM. The minimal form of the PQ-NMSSM, 
in which the singlino mass comes only from the singlino-Higgsino mixing and thus the LSP 
from the singlino-Higgsino sector becomes lighter than about 70 GeV, is shown to be tightly 
constrained by the Higgs invisible decay if the LSP is a candidate of the dark matter and 
it relic density is determined by the standard freeze-out process. 

Taking A = 0.7, which remains perturbative up to the PQ scale for tan /3 close to 1, the 
125 GeV Higgs can be obtained for stop mass around 500 GeV. Such a small tan (3 is favored 
as the Higgs invisible decay can be forbidden kinematically (?ti^o > m/i/2). However, it 
requires a large Higgsino component for the LSP and thus the recent XENONIOO bound on 
the nucleonic cross-section of the LSP at 90% C.L. excludes almost all the parameter region 
except a narrow band of the LSP mass close to the Higgs resonance point {fn^o — m/,/2 < 0.1 
GeV) which suppresses the thermal LSP relic density significantly. We also note that a 
late decay of the saxion inherent in the model can produce a huge amount of entropy 
diluting away the thermal LSP abundance. Then, the severe constraint on the LSP mass 
is invalidated. For larger tan /3, the LSP mass becomes smaller opening the Higgs invisible 
decay. In this case, a cancellation can be arranged to suppress the Higgs-LSP-LSP coupling. 
Combined with the charged Higgs mass bound at the LHC, the Higgs invisible decay 
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branching fraction is shown to become smaller than 10% for tan /3 > 9. In this parameter 
region, the NMSSM contribution to the Higgs boson mass is negligible and thus a large 
stop mass ~ 1 TeV is needed to get the Higgs mass of 125 GeV. As the LSP is almost 
purely singlino for large tan /3, its thermal freeze-out density is orders of magnitude larger 
than required. However, it turns out that a late decay of the saxion, washing out the 
dangerous thermal LSP relics again, can produce the right amount of non-thermal dark 
matter population for the singlino mass of order 10 MeV with A ~ 0.1. 

Generic collider signatures of the PQ-NMSSM with a light singlino-like LSP are multi- 
jets from multi-top/bottom and h/Z/W plus missing energy in the final states. Taking two 
benchmark points for small and large tan (3, we analyzed such signals at the 14 TeV LHC 
to find that a Sex discovery is possible for a few fb~^ of the integrated luminosity. 

We finally note that the unsatisfactory dark matter properties applied to the minimal 
PQ-NMSSM in the standard cosmology can be evaded in a more general PQ-NMSSM 
allowing a suitable bare singlino mass term. 
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Appendix 



A Neutralino mixing matrices 



The neutralino mass matrix is given in (HOI) , which can be diagonahzed by the method of 
perturbative diagonahzation as in 



(77) 
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(79) 



We spht the mass matrix into the diagonal part and the off-diagonal part, Ai = Ai^ + Ai^ , 
where = diag(Mi, M2, f^es, —f^^es, 0) and A^*^ is the rest. In the leading order, we have 



^(1) 



M 



o 



M 



D 



'■mm 



^0) 



In the perturbative limit, we find that 
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In the second order, 
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Then, we have 
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Keeping all these terms, we can write 
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For large tan/3, or more precisely, for Avtan/3//iofr > 1, A''i3 is very suppressed up to the 
second order. Thus, we include the third or higher order. The third order relation is given 
as 
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(3) (3) 

We here show only V53 and V54 since the other third order terms are neghgible. In the 
third order, we have 

V,f = ^f|^ + |flcos2/3(sin/3 + cos/3) 

■ sin 2/3 (sin /3 -cos (97) 



3 

eff 



■sin2/3(sin/3 + cos/3), (98) 



cff 



and 

Up to now, we have not sorted the eigenstates by magnitudes of eigenvalues. In the PQ- 
NMSSM, the singlino-like state is the lightest one in most parameter space. Therefore, 
we should change the mixing matrix in order to fit with the SUSY Les Houches Accord 2 
(SLHA2) convention i.e. N^i N^. 



B Higgs-bottom quark and Higgs-gauge boson cou- 
plings 

The CP-even Higgs mass matrix 0331) can be diagonalized by the mixing matrix S according 
to the SLHA2 convention. 

= S^jHJ'''^, (101) 

where H^'^^^ = [HdR, Hur, Sr) and H^^^^ are ordered in increasing mass. Then the various 
Higgs couplings can be expressed by the components of the mixing matrix. 

H.bLb'k : ^^.1, (102) 

Ha^^Z,: ^^,!M±M)(5^^cos/3 + ^,2sin/3), (103) 
2 

HiW^W- : g^,^{Sacosf3 + Si2smf3). (104) 
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C Higgs-neutralino couplings 



Referring to Appendix A in [3], the Higgs-neutrahno couphngs in the PQ-NMSSM are given 
by 

(105) 

where 

Each term can be explained in piecewise. The first term denotes ath Higgs couphng to 
ith and jth neutrahnos through the H^r component of the ath Higgs with the and S 
components of the neutrahnos, i.e. XH^rHu^ / \/2. The second and third terms are from 
permutations of the first term, i.e. XHurH^S / \/2 and XSrH^H^/ \/2. The fourth and fifth 
terms are from the SUSY U(l)y gauge couphngs of Hd and Hu muhiplets. Likewise, the 
last two terms are from the SUSY SU(2)l gauge interactions. 
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